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We study a relativistic spin-1/2 fermion subjected to a Dirac oscillator coupling and a constant 
magnetic field. An interplay between opposed chirality interactions culminates in the appearance 
of a relativistic quantum phase transition, which can be fully characterized. We obtain analytical 
expressions for the energy gap, order parameter, and canonical quantum fluctuations across the crit- 
ical point. Moreover, we also discuss the effect of this phase transition on the statistics of the chiral 
bosonic ensemble, where its super- or sub-Poissonian nature can be controled by means of external 
parameters. Finally, we study the entanglement properties between the degrees of freedom in the 
relativistic ground state, where an interesting transition between a bi-separable and a genuinely 
tripartite entangled state occurs. 
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I. INTRODUCTION 



the following Dirac equation 



Phase transitions describe an abrupt change in the 
physical properties of a system caused by the modifica- 
tion of its temperature. This phenomenon usually entails 
a change in the symmetry of the phases involved, and is 
commonly driven by thermal fluctuations. Consequently, 
classical phase transition cannot persist at zero tempera- 
ture where the absence of thermal fluctuations forbids the 
sudden change of phase. Nonetheless, other kind of fluc- 
tuations exist at zero temperature, the so called quantum 
fluctuations, which can also be responsible for a dramatic 
change in the properties of the system. In this case, the 
change is driven by the modification of certain couplings 
that describe the interactions between the microscopic 
constituents of the system, and is usually known as a 
quantum phase transition 

These critical phenomena arise in the thermodynami- 
cal limit of certain many-body systems, where the num- 
ber of particles TV — > oo. Usually, the description of such 
systems is extremely complex and one must deal with 
numerical methods. Nevertheless, there exist certain sit- 
uations where a simplified model, which can be exactly 
solved, captures the full essence of the problem display- 
ing such an abrupt change of the system properties. For 
instance, a collection of N two-level atoms interacting 
with a single mode of the radiation field, known as the 
Dicke model, displays a quantum phase transition whose 
features can be justified by means of a simple two-mode 
Hamiltonian in the thermodynamical limit 0, [1] . Other 
two-mode Hamiltonians which also display critical phe- 
nomena have been studied in the field of Quantum Op- 
tics [1, m , or in Nuclear Physics [6| . 

In this work, we shall be concerned with a relativis- 
tic toy model that involves two phonon modes and also 
displays a quantum phase transition. This critical phe- 
nomenon occurs in a fermionic relativistic harmonic os- 
cillator, also known as a Dirac oscillator 0, [1], when an 
additional constant magnetic field is applied. This rela- 
tivistic fermion of mass m and charge — e is described by 



ih 



dt 



= [ca (p - im/3wr + eA) -I- /3mc^] I*), (1) 



where j^f) stands for the Dirac 4-component spinor, r 
and p represent the position and momentum operators, 
w is the Dirac oscillator frequency, c stands for the speed 
of hght and f3 := diag(I, — I), :— off-diag((Tj, cTj) are 
the Dirac matrices related to the usual Pauli matrices 
Q. Here, the magnetic field is introduced by minimal 
coupling p — > p -I- eA , where A is the vector potential 
related to the magnetic field through B = VA A. On the 
other hand, the Dirac oscillator coupling is introduced by 
a non- minimal coupling p — im/Swr, where uj stands for 
the Dirac oscillator frequency. 

Here, we shall focus on a two-dimensional setup where 
the Dirac matrices become the well-known Pauli matri- 
ces: = crx,C(y = (Ty,/3 = ct^. In this scenario, Eq. ([T]) 
can be expressed as follows 



(Tj {pj ~ imPuzXj + eAj) + (Jzmc^ 



where is a 2-component spinor which mixes spin-up 
and -down components with positive and negative ener- 
gies. Remarkably enough, the Dirac oscillator coupling 
endows the particle with an intrinsic left-handed chiral- 
ity which is only present in this two-dimensional sce- 
nario [l3| . Conversely, the magnetic field coupling favors 
a right-handed chirality [llj , and therefore an intriguing 
interplay is set up. This system can be considered as a 
relativistic extension of chiral harmonic oscillators, which 
carry dual aspects of a certain symmetry (i.e. chirality 
in the plane). Such non- relativistic systems have been 
studied from a fundamental point of view [T^ , and an in- 
teresting connection to topological Chern-Simons gauge 
theories has also been pointed out [l^, . 

In this paper, we show how this relativistic chiral 
oscillator presents several intriguing critical properties, 
and offer an ideal scenario where to study the effect of 
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opposed chirality interactions in a two-dimensional set- 
ting. In Sec. [m we show how the relativistic Hamil- 
tonian in Eq. ([2]) can be exactly mapped onto a pair 
of simultaneous Jaynes-Cummings(JC) and Anti-Jayncs- 
Cummings(AJC) interactions jlSl . with right- and 
left-handed chirality respectively. This result differs sub- 
stantially from previous situations where only a distinc- 
tive chiral interaction appeared [lol . [lll | , and opens up the 
possibility to study a unique interplay between left- and 
right-handed interactions. In Sec. Illli we describe two 
limiting cases: a weak magnetic field regime where the 
system displays a remarkable left-handed chirality, and 
a strong magnetic field regime with an opposite right- 
handed chirality. For intermediate couplings, we obtain 
the complete energy spectrum and the associated eigen- 
states in Sec. IIVI To accomplish such task, we perform a 
unitary transformation which turns the Hamiltonian into 
a single-mode interaction which can be easily diagonal- 
ized. This exact solution shows how an unusual competi- 
tion between chiral terms arises, and leads to a relativis- 
tic level crossing phenomenon for a critical value of the 
JC and AJC couplings, which is described in detail in 
Scc.|Vl In this sense, several analogies with second order 
quantum phase transitions occur, such as the energy gap 
suppression, the divergence of quantum fluctuations, and 
the existence of an order parameter that reveals the ex- 
istence of a quantum phase transition. Furthermore, we 
also show that the statistical nature of the chiral phonon 
distribution displays an abrupt change across the critical 
point, where super-Poissonian chiral phonons turn into 
sub-Poissonian phonons and vice versa. Finally, in Ap- 
pendix |3 we give some details for the construction of the 
system eigenstates, which are closely related to SU{1, 1) 
coherent states. 



II. EXACT MAPPING ONTO A 
SIMULTANEOUS JC-AJC HAMILTONIAN 



the oscillator's ground state width associated to each fre- 
quency, and we have introduced i = x^yto account for 
the two possible directions of motion The substi- 

tution of these operators in Eq. followed by the in- 
troduction of the chiral annihilation operators for each 
frequency 



Or := -^{^x - ioy), ai 
a-r 75 (°^x - ioy), a/ 



Tflorc + ia,y), 



(4) 



and the consequent creation operators a}! , a| , aj! , aj! , leads 
to the following bichromatic relativistic Hamiltonian 



H 



mc 



1 



2^a 



-iV2ea/ + iy 2^4 -1 

^(5) 

where the two important parameters ^ Tiio/rac^, and 
^ :~ huj/mc^ represent the strength of the oscillator and 
magnetic field coupling with respect to the particle rest 
mass energy, respectively. This Hamiltonian in Eq. ([5]) 
can be rewritten as a simultaneous JC and AJC interac- 
tions with opposite chiralities 



H = 6a, 



(6) 



where S := mc^ stands for the detuning parameter pro- 
portional to the rest mass energy, H'^Q{gr) represents a 
right-handed Jaynes-Cummings Hamiltonian 



(7) 



with g,. :~ vmc? \j 2^ as the interaction coupling strength 

Analogously, the term H^j^igi) stands for a left-handed 
Anti- Jaynes-Cummings interaction 



gi(T+a] + g'la ai, 



(8) 



Let us first provide the exact mapping onto a pair of 
simultaneous JC and AJC couplings with opposite chiral- 
ities. We shall work in the axial gauge, where a constant 
magnetic field B = Be^ is described by the following vec- 
tor potential A := ^^ ^^^^ setup, the dynam- 
ics of a relativistic fermion is described by two different 
frequencies, the Dirac oscillator frequency uj and the cy- 
clotron frequency ujc ■— eB/m. In this regard, we must 
introduce a pair of sets of creation-annihilation operators 



with a similar coupling strength gi := imc\/2^. The 
relativistic Hamiltonian in Eq. ([6]) is depicted in fig. [U 
where the JC and AJC couplings give rise to a pair of 
channels through which the relativistic particle can per- 
form spin-flip transitions. 



III. WEAK AND STRONG MAGNETIC FIELD 
REGIMES 



1 fl .A \ 



1 / 1 



.A 
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1 / 1 



.A 



/f a"' 'J''' ' 



(3) 



associated to the frequencies to and lj ;= loc/2 respec- 
tively. Here, A := y^h/muj and A :— ^ h/mCo represent 



The exact mapping of this relativistic model onto a si- 
multaneous JC and AJC couplings allows a neat descrip- 
tion of the interaction in terms of a energy level scheme 
(see fig. [1]). The presence of the pair of simultaneous 
couphngs H^j^{gi), and Hf(^{gr) in Eq. ([5]), forbids a 
straightforward solution of the complete Hamiltonian by 
analogy to JC-like solutions [l3j[II|- Nonetheless, these 
techniques can be used to obtain the energy spectrum 
and corresponding eigenstates in two limiting regimes: 
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Weak magnetic field ^ <C ^ 



In this regime, the Hamiltonian in Eq. ([5|) becomes 



= 25a, + gia+a] + .gfa-fl; + OU/i], (9) 



and we can obtain a first approximation to the energy 
spectrum neglecting corrections of order O {(,/£,). Con- 
sequently, the relativistic Hamiltonian can be approxi- 
mately described by a left-handed AJC coupling which 
can be directly solved. The energy spectrum becomes 
^ ^ ±mc^[i + 2^{ni + l)Y^^, where = 0,1... 

represents the number of left-handed phonons. The cor- 
responding eigenstates form an AJC-doublet which can 
be described as follows 

\±E:^)= CZ,±\ni + 1)|xt) TiC„-„^K)IXi), (10) 

where we have introduced the set of left-handed Fock 
states \ni) := (n/!)~-'^/^(aJ)"' |vac) , and certain normal- 
ization constants C'^^ j^ := [{El";^ ± mc^)/2Ei;f/^ . 



B. Strong magnetic field (, ^ ^ 

In this limit, the Hamiltonian in Eq. ([5]) becomes 

= 25a, - gra+hr - g^-^ + q (c/c) , (H) 

and we can find the energy spectrum up to 0(^7^), which 
becomes E" = ±EI^^ = ±mc^[l + 2^{nr + 1)]^''^, where 
fir = 0, 1... is the number of right-handed phonons. The 
corresponding eigenstates form a JC-doublet with the fol- 
lowing structure 

I ± El) = Ct,±|n,)|xT) T ^IJnr + l)|xi), (12) 



lXT> 




FIG. 1: Level scheme representation of relativistic interac- 
tion, which gives rise to spin-flip transitions accompanied by 
the simultaneous creation and annihilation of chiral phonons. 
Note that there are two different spin-flip channels depending 
on the phonon chirality . 



with \hr) := (n-rO "'^/^(aj^)"'' |vac) as right-handed Fock 
states, and C|^ ± := [{Ei^ ± mc'^)/2E?^f^'^ . 

The latter results can be understood as follows, in the 
regime where ^/^ — > 0, the relativistic system is found in 
a phase with left chirality. Conversely, in the limit where 
^/^ cx), a notorious right-handed chiral phase arises. 
Therefore, a modification of the coupling strengths ^, ^ 
leads to a change of the chiral symmetry, which can be 
interpreted in the language of quantum phase transitions. 
We will consider a zero-temperature setup, where ther- 
mal fluctuations do not exist and the chirality transition 
can only be driven by quantum fluctuations. As we show 
below, the complete relativistic energy spectrum can be 
obtained for all possible couplings ^, ^. This spectrum 
presents two crucial properties which resemble the usual 
setting in quantum phase transitions: the energy spec- 
trum is non- analytical for a critical coupling (C/Oc = 1; 
and the system becomes gapless at this critical point 
A£; ^0 when 1(1/0 -d/Oc I ->0. 



IV. EXACT SOLUTION: ENERGY SPECTRUM 
AND ASSOCIATED EIGENSTATES 



The particular structure of the Hamiltonian in Eq. ([5|) 
suggests a different description where the spinorial de- 
grees of freedom are coupled to a unique bosonic opera- 
tor b(x{ai- (C/0^^^4) which mixes both chiralities de- 
pending on the relative coupling strength ^/^. In this re- 
gard, wc derive a unitary transformation which converts 
the bichromatic Hamiltonian in Eq. ([6|) into a monochro- 
matic JC (AJC) term that involves such bosonic degree of 
freedom b with a certain chirality that depends on exter- 
nal parameters ^ > ^ (^ < ^ . This transformation, which 
turns out to be a two-mode squeezing operator, allows an 
insightful derivation of the energy spectrum and its as- 
sociated eigenstates whose properties depend strongly on 
the magnitude of the magnetic field applied. As we shall 
discuss, certain properties of the system, such as chiral- 
ity, squeezing, phonon statistics, and entanglement, are 
conditioned to the value of the magnetic field. Moreover, 
as the magnetic field is varied, a sudden quantum phase 
transition occurs at ^ = ^ where both chiralities con- 
tribute identically, and the critical theory becomes that 
of a relativistic free fcrmion. 



A. Left-handed regime ^ < ^ 

Under these circumstances, the relativistic Hamil- 
tonian in Eqs. ([6])- ([8]) can be unitarily mapped onto 
a single-mode Anti-Jayncs-Cummings Hamiltonian by 
means of the following unitary transformation 



(13) 
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where the real parameter a depends on the relative 
strength of the oscillator and magnetic field couplings 



a:.iarctanh(^^ 



(14) 



and we have introduced /i := (A/A — A/A)/2, and 
A := + The transformation laws for the chiral 



operators under the unitary in Eq. (|T3|) can be described 
as follows 

Uj,a/UQ = ^cosh(AQ:) + ^sinh(Aa)^ a; — — sinh(Aa)a]^, 

UljOrUa = ^cosh(Aa) — ^sinh(Aa)^ a,. — — sinh(Aa)aJ, 

(15) 

which lead to the following transformation of the com- 
plete original Hamiltonian into a single-mode AJC term 

H^cial) V^HVi = 5a, + g\n+a\ + [g'lYn-au (16) 

with a modified coupling strength g[ :~ imc^-v/S^, where 
Ci{i,Cj := ^ - f - 2^(1^)1/2 is related to the initial rel- 
evant parameters. Once the single-mode Hamiltonian in 
Eq. (fTH]) has been obtained, we can calculate the energy 
spectrum and corresponding eigenstates of the original 
Hamiltonian by solving the single mode coupling. Two 
Hermitian operators related through a unitary transfor- 
mation share a common spectrum, and therefore the en- 
ergy eigenvalues of the Hamiltonian in Eq. ([5]) can be 
obtained in analogy to the weak magnetic field regime 



(17) 



where n; = 0,1... represents the number of left-handed 
quanta. The associated eigenstates arc obtained by 
means of the unitary transform applied to the single- 
mode left-handed eigenstates 

|±i?„,) = Ul {Cn,,±\ni + 1)|xt) T iC„,Tl"/>IXi» - (18) 

where the normalization constants depend on the ener- 
gies C„,^± := [(i?„, ±TOC^)/2£'„,]^/^. The transformation 
in Eq. (|13p can be rewritten in the monochromatic sce- 
nario as 



(19) 



where jl := (A/A-l- A/A)/2. This transformation can be 
immediately related to a two-mode squeezing operator in 
the context of Quantum Optics with squeezing parameter 
z := — a/i/2 € M [l^. The action of such an squeezing 
operator over left-handed chiral Pock states gives rise 
to SU(1,1) coherent states \z,ni) U„|n;) |vac)r [s^fsee 
Appendix 1^ for some details) 



\z,ni)^Afn, £ \ I ^ ' (~ l)"tanh'" \z\\m + ni,m), 



m—O 



ni\m\ 



where we have introduced the normalization constant 
Mm ■= cosh~^"'^^''|z|. Since right-handed operators are 
not present in the effective Hamiltonian of Eq. (fT6|) , and 
consequently do not participate in the relativistic dy- 
namics, we have chosen the right-handed vacuum for 
simplicity. The state corresponding to ri; = can be 
identified with the two-mode squeezed vacuum state, an 
archetypical state in the field of Quantum Optics that 
becomes highly non-classical for a large squeezing param- 
eter. Therefore, the energy eigenstates of the relativistic 
fermion in Eq. (|18p can be expressed as 

|±S„,) =C„„±|z,n, + l)|xT)TiC„,,^|z,n,)|xi), (21) 

which is identical to an AJC-like doublet where SU(1,1) 
coherent states are entangled with the relativistic spinors. 
Such an state possesses remarkable non-classical proper- 
ties, such as spin-orbit entanglement or sub-Poissonian 
statistics. With respect to the weak magnetic field eigen- 
states in Eqs. ([T0|) . these states present certain novel 
features, such as an inter-mode chiral entanglement, or 
bosonic statistics which depend on the coupling strength 
^/^. All these interesting properties will be described in 
detail in forthcoming sections. We shall be also interested 
in the properties of the fermionic ground state Eg = mc^ 



\9) = k,0)|xT) = -^^f2{-iy''te^nh^\z\\m,m)\x^), 
cosh z '^-^ 

' ' m=0 

(22) 

which can be notably interpreted as a spin-up squeezed 
vacuum state, where the squeezing parameter z = — Q;/i/2 
depends on the relative coupling strengths 



B. Right-handed regime ^ > ^ 

In this case, the original Hamiltonian in Eqs. ([H])-® is 
transformed into a Jaynes-Cummings Hamiltonian under 
the action of 

Us :=e"('"°--">"'), (23) 
where the parameter a in this regime becomes 



a 



arctanh 



(24) 



(20) 



The chiral operators are once again transformed accord- 
ing to Eqs. psp when the substitution a — s- a is per- 
formed. In this situation, the transformed Hamiltonian 
becomes a single-mode Jaynes-Cummings term 

H^cid'r) Uf^i/Ut = da, + g',a+ar + ig'r)*a-al, (25) 

where the new coupling strength is g'^ := mi(?\f2^^ and 
Cr(^,0 := C - ? + 2m(C0^^^- Analogously to the left- 
handed regime ^ < we obtain a single-mode Hamilto- 
nian which can be easily diagonalized following the same 



5 



procedure as in the strong magnetic field limit, and pro- 
vides the solution to the complete bichromatic interaction 
in Eq. ([6]). In this sense, the energy spectrum becomes 

E = ±En^ = ±mc^ ^J^ + 2CritOinr + l), (26) 

where = 0, 1... represents the number of right-handed 
quanta. Applying the unitary transform in Eq. (|23p to 
the single-mode eigenstates, one obtains the correspond- 
ing cigcnstates of the complete original Hamiltonian 

I ± En^) = ^1 {Cn,^,±\nr)\xi) T ia^M^r + l)\xi)) , 

(27) 

where Cn^^± := [{En^ ± mc'^)/2En^^]^/^ . In this regime, 
the transformation in Eq. (|23p in the single-frequency 
domain becomes 

=e-^(-°'-"'+°^"'), (28) 

which can be reinterpreted once more as a two-mode 
squeezing operator with a different squeezing parame- 
ter z := —afl/2 € R. This operator transforms the 
right-handed chiral Fock states \hr) into SU(1,1) coher- 
ent states |i, nr) := Ut |vac);|nr), where we have chosen 
the left-handed vacuum for simplicity (see Appendix |^ 



\~z,nr) ^Mn^ V / + ) ! („ ^ )r»tanl/" I g 1 1 m+n. , m) , 

m=0 V 

(29) 

where we have introduced the following normalization 
constant Af^^ :— cosh~'-"'''''^''|i|. Such coherent states in 
Eq. (|^ appear in doublets in the fermionic eigenstates 
of Eq. (EH) 

|±-Es„) ^Cf,^,±\S,fir)\x])TiCn,.,^\S,hr + l)\xi)- (30) 

Once again, one can identify these eigenstates as JC-likc 
doublets where SU(1,1) coherent states become entan- 
gled with the relativistic spin degree of freedom. Con- 
versely to the fermionic ground state of the left-handed 
regime in Eq. ([22]) . the ground state in this regime 
Eg = mc^Vl + 2Cr, cannot be expressed solely by the 
two-mode squeezed vacuum, but rather by 

\g)^Co+\~z,0)\x^)^iCo-\~z,l)\xi), (31) 

which besides the two-mode squeezed vacuum, also in- 
cludes the fir = 1 SU(1,1) coherent state. 

C. Critical regime ^ = ^ 

In this regime, the effective Hamiltonian can be di- 
rectly obtained from Eq. ([2]) 

2 

iJ2^^5<= = c ^ Gjpj + a,mc2 , (32) 



which corresponds to the Dirac Hamiltonian of a two- 
dimensional free fermion. In the critical regime ^ = ^, 
the magnetic field coupling cancels the effect of the Dirac 
string coupling, and the fermion behaves as a free rela- 
tivistic particles. The critical energy spectrum becomes 

E'' = ±El = ± V"i2c4 + p2c2, (33) 

where p = {pa:,Py) stands for the two-dimensional 
fermion momentum. The corresponding eigenstates can 
be described as follows 

l mc^±E- / c(p.+ip,) A 

' ^ - V ±2i?c (,I^T> + ^,2^eJxi) ) IP>. 

(34) 

where |p) := \Px,Py) are two-dimensional plane wave so- 
lutions. Note that these solutions describe the relativistic 
fermion at the critical point ^ = 5, but any small pertur- 
bation of the magnetic field or the Dirac string coupling 
will dramatically change the system properties. As we 
discuss below, this dramatic change shares many analo- 
gies with a quantum phase transition. 

We have thus provided a complete solution of the rel- 
ativistic Hamiltonian that describes the properties of a 
Dirac oscillator subjected to an additional constant mag- 
netic field. A unitary transformation that connects the 
bichromatic full Hamiltonian with single-mode JC-like 
interactions has been described in Eqs. (|23p . With 

the aid of such transformation, one can obtain the exact 
energy spectrum for all possible values of the coupling 
parameters ^, ^, which responds to the analytical expres- 
sions in Eqs. PT)) . (P5|) . In this regard, one can study 
the transition between the weak ^/^ — > and strong 
^/^ — > oo coupling regimes which endows the relativistic 
system with a chiral symmetry change, as discussed in 
the previous section. 

From a practical point of view, the magnetic field de- 
pendence of this relativistic system might be understood 
as an accessible method to prepare a certain relativistic 
state with specific properties. An adiabatic control of the 
applied magnetic field opens up the possibility of control- 
ling a wide range of properties of this relativistic fermion, 
such as its chirality, squeezing, phonon statistics, and en- 
tanglement, which shall be discussed below. 

V. CHIRALITY QUANTUM PHASE 
TRANSITION 

In previous sections, we have described how the value 
of the relative coupling strength ^ is responsible for the 
chirality of the system, which can be right- or left-handed 
by an appropriate tuning of the coupling strengths. This 
drastic modification can only be driven by means of quan- 
tum fluctuations, and therefore a quantum phase transi- 
tions occurs in this relativistic system. In this section we 
study several properties which are clear signatures of a 
quantum phase transition phenomenon, such as the van- 
ishing of the gap at the critical point, the existence of 
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an order parameter which takes on different values in the 
distinct chiral phases, the divergence of quantum fluctu- 
ations and the maximization of entanglement across the 
critical point. 



A. Energy level crossing 

Here, we study the properties of the relativistic energy 
spectrum described in Eqs. (|17p . (pS)) . which are repre- 
sented in fig. [2l Here we represent the energies of dif- 
ferent eigenstates | ± -En,), | ± En^) with respect to the 
relative strength ^/^. 



We clearly observe two crucial signatures of a quantum 
phase transition (see also fig. [3] for a two-dimensional sec- 
tor of the above energy spectrum): 

• The energy spectrum is non-analytical for a critical 
coupling = 1, 

• The system becomes gapless at this critical point 
Ai? ^0 when 1(1/0 -d/Ocl^O. 

Furthermore, we can obtain an analytical expression 
of the gap close to the critical point. In terms of the 
couphng strengths gr and gi in Eqs. ([7])- ([8]), we obtain the 
following universal scaling law for the energy gap close to 

~ 1/2 

the critical point {gr/gi)c = {£./£,)c = 1 




5r _ 






9r 
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m 






m 


\9lJc 



and we can readily identify the scaling exponents zv ^ I. 
The critical exponent z is a dynamical exponent related 
to the decay of the system fluctuations with time, and in 
a Lorentz invariant model is expected to be 2 = 1. Con- 
sequently, the remaining exponent, which characterizes 
the scaling of the correlation length close to the phase 
transition, becomes = 1 in this relativistic model. 



FIG. 2: Fermionic energy spectrum as a function of the rela- 
tive coupling strengths ^/^ and the number of chiral phonons 
ria ijia = n; if ^ < ^, and ric, = fir ^ > £,)■ The two energy 
sheets correspond to positive- and negative-energy solutions 
I ± Er,,) for I < e, and | ± Ef,^) for f > 



E/mc^, 




FIG. 3: Energy levels for the first excited eigenstates as a 
function of the relative coupling strengths showing the 
non-analytic behavior at (5/C)c = 1. 



1.4 - 
1.2 - 



FIG. 4: Ground state energy as a function of the relative 
coupling ^/^. As the magnetic filed is increased, the ground 
state energy changes from Eg = m(? Eg — mc^y^l + 2C,r, 
which becomes non-analytical at the critical point. 



Finally, it is also important in quantum phase transi- 
tions to study the properties of the system ground state. 
In this case, the fermionic ground state in the left-handed 
regime is described by Eq. ((22|) . whilst the ground state 
in the right-handed scenario becomes that of Eq. ([31]) . 
The energy of such ground state is non-analytical at the 
critical point a follows from fig. 01 
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B. Order parameter 

Another clear signature of quantum phase transitions 
is the appearance of an order parameter, a physical mag- 
nitude which acquires different values in the phases in- 
volved and becomes indeterminate at the critical point. 
Therefore, an order parameter is able to witness the 
abrupt change in the properties of the system. In this 
relativistic scenario, the order parameter turns out to 
be the z-component of the orbital angular momentum 
Lz = xpy — ypx, an observable which is intimately re- 
lated to the symmetry properties of the system, since 
Jz = Lz + Sz is a conserved quantity. Note that this 
operator is also connected to the system chirality, since 
left-handed states fulfill {Lz)i < whereas right-handed 
states fulfill {Lz}r > 0. In order to evaluate the angu- 
lar momentum expectation value, it is useful to express 
the unitary transformations in Eqs. (flQl) . ([28]) in terms of 
canonical conjugate position and momentum operators 



(36) 



Using the transformations in Eqs. one can obtain 

the expressions of the transformed position and momen- 
tum operators 



UqXjUJj = e 2 a;i, UaPiUj, = e pi, if £.<£,, 



2 Pi, if e > 



(37) 



with i ~ x,y, which leads to the corresponding relations 
for the orbital angular momentum 



UaLzUj, =Ua(a;py-?/p^)Uj, =Lz, if ^ < 
UfiLzUt =Ufi(xp,y-2/p,)Ut =Lz, if |>e 



(38) 




(C/C)c = 1 



FIG. 5: Mean value of the z-component of the orbital angular 
momentum across the critical point. Note that, as the critical 
region is traversed, the sign of this order parameter changes 
and indicates that the phase transition has occurred. 



We observe from the expressions in Eqs. (|38|) that the 
z-component of the orbital angular momentum operator 



is not altered by the squeezing transformation. In the 
language of Lie algebras, this property is equivalent to 
Lz being the Casimir operator associated to the underly- 
ing SU(1,1) algebra (see Apendix[X|. Therefore, we can 
easily obtain its expectation value in the eigenstates of 
Eqs. ((I8]),(l27l) 



"""'^""'"l^^O, if |<^, 



{Lz)\±E^^) = -h{ni + 



2E„ 



(L 



= +n I + ) > 0, if |>C, 



(39) 



which show how the orbital angular momentum takes on 
negative values in the left-handed chiral phase, whereas it 
attains positive values in the right-handed chiral regime 
(see fig.[5|). Consequently, the orbital angular momentum 
plays the role of an order parameter which witnesses the 
quantum phase transition and macroscopically reveals 
such effect by a change of its sign. 



C. Divergence of quantum fluctuations 

In the vicinity of a critical point, the gap of the sys- 
tem becomes negligible and excitations can be easily pro- 
duced (see fig. Under such conditions, the system 
becomes highly fluctuating. As discussed previously, a 
quantum phase transition can only be driven by quan- 
tum fluctuations which lead to an abrupt change in the 
physical properties of the system when the critical region 
is crossed. In our case, it is possible to calculate analyti- 
cally the divergences in the quantum fluctuations of the 
fermion position Axi ~ \J {xf) — (xi)'^. 

We shall not only focus in the ground state fluctu- 
ations in Eqs. (P^ . I^Tj) . but wc shall investigate the 
quantum fluctuations of the whole energy eigenstates in 
Eqs. (PT|) . ([501) . Using the transformations described in 
Eq. ((371), we obtain 



Axi 



Axi 



|±B„,) 



Axi 
Ax, 




(40) 



where i = x,y, Axi\vac '■= A/^/2 stand for the fluctua- 
tions of the vacuum in modes ar,ai, whereas Ax^jvac •= 
A/V2 arc the vacuum fluctuations in modes dr,ai. We 
have also introduced r]± :— ni + 3/2 ± mc^/2£'„, , and 
?7-i- :— fir + 3/2 ± mc^ /2En^, which depend on the num- 
ber of quanta. We have represented in fig. [6] the fluctu- 
ations in the fermion position for different eigenstates as 
the relative coupling ^/^ is varied. We observe how these 
fluctuations diverge at the critical point (C/Oc = 1 for 
all the eigenstates, representing yet another signature of 
a quantum phase transition. 

Nonetheless, we have remarked in Eqs. p^ . (|28p that 
the unitary transformation involved corresponds to a 
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FIG. 6: Quantum fluctuations in the fermion position l\Xi for 
the first excited positive-energy eigenstates | +Eni ) , \ +En^ ) as 
a function of the relative coupling strengths ^/^. We note that 
the negative-energy eigenstates display similar fluctuations in 
position. 




FIG. 7: Quantum fluctuations in the fermion momentumAp, 
for the first excited positive-energy eigenstates | + _B„, ) . | + 
Ef,^} as a function of the relative coupling strengths ^/^. We 
note that the negative-energy eigenstates display similar fiuc- 
tuations in momentum. 



squeezing transformation in the language of quantum op- 
tics. As a consequence, a squeezing of certain fluctuations 
must also become apparent for a certain observable. This 
is the case of the fermionic momentum, whose fluctua- 
tions Api = y/ (pf) — {pi)"^ for different eigenstates are 
the following 



Ap^\^E„^^ = Ap^\^^^y/^;:;Je if c<?, 

^ Api if ^>i, 



Ap, 



(41) 



where Apjjvac = h/ A\/2^ and Ap^jvac = h/Ay/2 stand 
for the vacuum fluctuations in each regime. Remarkably, 
the fluctuations in the momentum vanish as the system 
approaches the critical point (sec fig. E]). Therefore, the 
squeezing occurs in the rclativistic momentum. 

With these observations, we have a full description of 
the system for any value of the relative coupling: 

• Left-handed regime ^ < The relativistic 
eigenstates arc described by an AJC-like doublet 
of SU(1,1) coherent states entangled with Pauli 
spinors (|2ip . which are rotating clockwise (see 
Eq. ^). 

• Critical regime ^ = The critical properties 
are described by those of a relativistic free fermion. 

• Right-handed regime ^ > The relativis- 
tic eigenstates arc described by an JC-like dou- 
blet of SU(1,1) coherent states entangled with Pauli 
spinors (j30[) . which are rotating counterclockwise 
(see Eq. (P])). 



D. Phonon statistical properties 

In this section we study the statistical properties of 
the chiral phonon distribution. The statistical properties 
of a bosonic ensemble can be classified according to its 
quantum fluctuations: 

• Poissonian statistics.- In this case, the phonon 
number distribution is a Poissonian random vari- 
able. Quantum fluctuations in the number of 

1^ a is the 



bosons fulfill An = \/{n), where n = a'a 
bosonic number operator, and An = (n^) — (n)^ 



• Super- Poissonian statistics.- In this regime, the 
quantum fluctuations of the bosonic ensemble are 
bigger with respect to the Poissonian distribution 
An > a/ (n), and therefore the phonons are said to 
be noisier. 

• Sub-Poissonian statistics.- The quantum na- 
ture of a bosonic ensemble allows a further pos- 
sibility, which is an evidence of non-classical be- 
haviour, that of sub-Poissonian statistics. Under 
these circumstances, the statistical fluctuations be- 
come lower than those of a Poissonian distribution 
An < \/(n) , and consequently the bosons are said 
to be quieter quantum entities. 

In this section we show how the relativistic phonon 
distribution attains different statistical properties de- 
pending on the relative coupling parameter ^/^. In 
particular, we show how the quantum noise of the chi- 
ral phonon distribution changes from sub-Poissonian to 
super-Poissonian as the critical point is crossed. In order 
to quantify the Poissonian character of the chiral boson 
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ensemble, we shall make use of the chiral Mandel Q pa- 
rameters, defined as follows 



1, 



h = —, — ^ 1- 



(42) 



A positive sign in such parameters reveals the super- 
Poissonian nature of the bosonic ensemble, whereas a 
negative sign shows the sub-Poissonian bosonic statis- 
tics which stress the non-classical nature of the ensem- 
ble. The corresponding fluctuations in the system eigen- 
states (PT|) . ([50]) ■ can be calculated using the expressions 
in Eqs. fXe )) . (|T7)) in Appendix In the left-handed 
regime ^ < ^, we have found 



{nr)\±E^^) = ?y±sinh |z|, 
("■/>|±£;„,) = ?7±cosli^|z| - 1, 
An^l^j.^ ^ = ?7±sinh^|z|cosh^|z| -|- Ksinh*|z|, 
Ani|^j_^^ ^ = r7±sinh^|z|cosh^|z| -f KC0sh^|z|, 



(43) 



where the parameters •r]± have been previously intro- 
duced, and K := — rv?c^ j E^^. Analogously, in the 
right-handed case ^ > ^ 



r/|±£;-J - f\^COS\?\l\ - 1, 

ryipSinh^|z|cosh^|z| -|- KCOsh^|z|, 

2|~| 1 2|~| , ~ • 1 4|~| 



{fii 



(44) 



Ani|^_i_jj_ ^ = T^zpsinh^ |z|cosh^|z| + Ksinh"'|z 



where f}± have already been introduced, and the parame- 
ter k := |(1 — Tn^c^ / E^^). These expressions allow us to 
study in detail the statistical nature of the chiral phonon 
ensemble by means of the Mandel parameters in Eq. 
which have been represented in figs. [8] and [9] In this fig- 
ures we notice the following: 

• Left-handed regime ^ < In this limit, the 
chiral Mandel parameters become: 

Qi < Q =^ Sub-Poissonian statistics of quieter left- 
handed phonons, 

Qr > =i>Super-Poissonian statistics of noisier 
right-handed phonons. 

• Right-handed regime ^ > In this limit, the 
chiral Mandel parameters become: 

Qi > =>Super-Poissonian statistics of noisier left- 
handed phonons, 

Qr < =J> Sub-Poissonian statistics of quieter right- 
handed phonons. 

Therefore, we observe how the statistical nature of the 
relativistic phonons is completely controlled by the rela- 
tive coupling strength. Remarkably enough, the level of 
quantum noise in the chiral ensembles can be controlled 




FIG. 8: Left-handed Mandel parameter Qi for the energy 
eigenstates {| ± En^), \ ± En^)} as a function of the relative 
coupling strengths ^/^. Note that as the magnetic field is 
raised, the ensemble of left-handed phonons changes from sub- 
to super-Poissonian statistics. 




FIG. 9: Right-handed Mandel parameter Qr for the energy 
eigenstates {| ± En^), \ ± En^)} as a function of the relative 
coupling strengths ^/^. Note that as the magnetic field is 
raised, the ensemble of right-handed phonons changes from 
super- to sub-Poissonian statistics. 



by means of external parameters. Furthermore, a notori- 
ous transition between sub- and super-Poissonian statis- 
tics occurs as the system crosses the critical region. This 
fact shows how a quantum phase transition can also be 
intrinsically related to the statistical nature of the sys- 
tem. Note also that close to the critical point Qi,Qr ^ 1, 
and both ensembles become super-Poissonian indicating 
how quantum fluctuations diverge at the critical point of 
a quantum phase transition. 
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E. Entanglement at the critical point 

In previous sections, we have seen how several signa- 
tures of a quantum phase transition occur in this rela- 
tivistic system. We have described in detail the closure 
of the energy gap at the critical point, the existence of 
an order parameter which turns out to be the orbital an- 
gular momentum, and the divergence of quantum fluctu- 
ations as one explores the critical region. An additional 
feature in quantum phase transitions is the divergence 
of quantum correlations as the system crosses the critical 
region Correlations of a quantum nature, known 

as entanglement, lie at the heart of Quantum Mechanics 
and arc responsible of non-local phenomena in the quan- 
tum regime. Furthermore, they are of outmost relevance 
in the fields of Quantum Information and Computation, 
where the y co nstitute a resource for information process- 
ing tasks I21I. [2^. [23} . In the particular case of fermions, 
entanglement has been previo usly studied for relativis- 
tic field theory [H, [2^ (see [Hi for a recent review). 
In this section we provide a description of fcrmionic en- 
tanglement across a quantum phase transition, that of a 
relativistic Dirac oscillator subjected to a magnetic field. 




FIG. 10: The relativistic degrees of freedom can be described 
as a tripartite system with a Hilbert space Ti — Hr ^TLi® 
Ha that includes continuous-variables corresponding to the 
chiral degrees of freedom, and discrete variables related to 
the spinorial degree of freedom. 

In our case, the underlying Hilbert space can be de- 
scribed as a tripartite system composed of continuous- 
variables associated to the chiral degrees of freedom, and 
discrete- variables related to the spin degree of freedom 
(see fig. nop . To quantify quantum correlations in this 
hybrid system we make use of the machinery of discrete- 
variable [131 and continuous-variable [l^ entanglement 
measures. In particular, wc obtain the von Neumann en- 
tropy of the reduced density matrices over all possible 
bipartitions 



Pi 

Pr 
Ps 



Tr,(Tr,|V')(^|), 
= Tr, {TTi\i;){^\), 
:Tr, (Tr,|V')(^|), 



(45) 



where the subindcxes /, r, s stand for the left-handed, 
right-handed, and spinorial degrees of freedom, \ip) is a 
particular pure state, and the von Neumann entropy is 



defined as S{p) := — Tr(plogp). Such magnitude, when 
calculated over reduced density matrices, quantifies the 
amount of entanglement between two parties. In our 
case, it quantifies the amount of entanglement between 
different relativistic degrees od freedom. 

Left-handed regime ^ < The ground state of 
the relativistic spin- 1/2 oscillator in Eq. is described 
by means of a spin-up two-mode squeezed vacuum state, 
with the following reduced density matrices 



Pi 



cosh 



^ 00 

-2— ^ tanh^"'|z||n/)(n;|, 



ni=0 

00 



Pr — ^ — tanh^""^ |z||?T.r)(nr 

cosh^lzl ' " 

' ' 7lr — 

Ps = IxtXxtI- 



(46) 



The reduced chiral degrees of freedom are described by 
means of a mixed thermal state with an effective temper- 
ature Teff := hLu/2kBlog(coth\z\), where ks is the Boltz- 
mann constant. Remarkably, the entropy of a thermal 
state can be given in terms of its mean number of quanta 
(n) [ii,!!^, and we obtain Si -.^ S{pi) ^ sinh2|z|log(l -f- 
cosech^|z|) -t- log(cosh'^|2|) = S{pr) =: Sr- There exists 
thus a certain amount of chiral entanglement, which de- 
pends on the relative coupling strength and diverges as 
one approaches the critical point (C/Oc = 1 ( sec fig. [TT|) . 
Conversely, the reduced spin density matrix is in a spin- 
up pure state with vanishing entropy Ss := S{ps) = ( 
see fig. [T^ . This remarks the fact that there is no spin- 
orbit entanglement between the discrete- and continuous- 
degrees of freedom, which is in accordance to the ground 
state structure in Eq. ([^^ . which is biseparable with re- 
spect to the spin-orbit bipartition. 



S(Pa 








— — ■5''- 



1 1.5 2 

FIG. 11: Von Neumann entropy of the single-mode reduced 
states Pa, with a — l,r indexing the mode chirality, as a func- 
tion of the coupling strength ratio ^/^. One directly observes 
that the entanglement of the continuous-variable degrees of 
freedom with the rest of the system diverges as one approaches 
the critical region. 
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Right-handed regime ^ > In this limit, the rel- 
ativistic ground state in Eq. ([51]) can no longer be de- 
scribed by a spin-up squeezed vacuum state, but rather 
by means of a SU(l,l)-doublet where chiral coherent 
states are entangled with spin states. The reduced den- 
sity matrices become 

PI = -^JTZ-. E tanh^'^'l^l 7+ + \ni){nil 
cosh \z\ V smh \z\J 

Pr = 727^ y] tanh^"'-|z| (7+ + ''"■^aT-i^^ ) \^r){nr\, 

cosh \z\ g V cosh \z\ J 

Ps ^ 7+Ixt)(xtI +7-IXi)(Xil, 

(47) 

where 7-1- := {y/T+'2Q^ zL l)/2vT+"^Cr- In this regime, 
the reduced chiral degrees of freedom cannot be described 
by means of effective thermal states, nor can the reduced 
spin degrees of freedom be identified with a pure state. 
Nonetheless, the density matrices in Eqs. (|47)l are already 
in diagonal form, and the von Neumann entropy can be 
expressed as follows 

00 

Si ^ - eiiogQ^, 

00 

= h (safe)) + ^log (^±f)" ' 

(48) 




FIG. 12: Von Neumann entropy of the spinorial reduced state 
Ps as a function of the coupling strength ratio ^/^. We observe 
that the entanglement of the spin degree of freedom to the rest 
of the system increases with the magnetic field. 



tudcs. which represent entanglement between all pos- 
sible bipartitions of the system have been represented 
in figs. [11] and [T^l We observe certain peculiarities 
with respect to the left-handed regime, such as a finite 
spinorial entanglement which tends to its maximal value 
^max _ jQg2 magnetic field is raised ^/^ ^ 00. 

We also observe that the continuous-variable entangle- 
ment is no longer equal for opposed chiral modes, which 
is a consequence of a finite orbital angular momentum in 
the ground state. In addition, one can see from Eq. (|3ip 
that the relativistic ground state presents genuine tripar- 
tite entanglement, where all the degrees of freedom are 
correlated among each other. 



VI. CONCLUSIONS 

In this work we have studied the properties of a rel- 
ativistic spin-1/2 Dirac oscillator subjected to a con- 
stant magnetic field. The relativistic Hamiltonian can 
be mapped onto a couple Jaynes-Cummings and Anti- 
Jaynes-Cummings terms, which describe the interaction 
between the relativistic spinor and bosons that carry an 
intrinsic chirality. These models, which are of paramount 
importance in Quantum Optics, become useful in this 
relativistic scenario, and allow an insightful description 
of interplay between opposed chirality interactions. The 
study of this relativistic system in the limits of a weak 
and strong magnetic field points towards the existence of 
two different phases, each described by an opposed chi- 
rality. In the intermediate regime, an intriguing trade-off 
between left- and right-handed interactions leads to the 
appearance of a quantum phase transition, which we have 
described in full detail. 

The remarkable possibility to find the exact solution of 
this relativistic system has allowed us to describe analyt- 
ically several properties of this phase transition. We have 
discussed the non-analytic nature of the energy spectrum 
at the critical point, and the sudden extinction of the en- 
ergy gap, which are archetypical properties of quantum 
phase transitions. We have also proved that there is an 
order parameter, the 2-component of the orbital angu- 
lar momentum, that witnesses the chiral phase change. 
In the critical region, we have shown that the quantum 
fluctuations in the fermion position diverge, which is a 
typical sign of quantum phase transition. Conversely, the 
quantum fluctuations in momentum vanish at the criti- 
cal point, since a particular squeezing of this canonical 
variable occurs due to the interplay of the interactions. 
We have discussed the bosonic ensemble statistics, and 
their behavior across the phase transition. In this re- 
gard, we have found that super- or sub-Poissonian nature 
of such ensemble can be controlled at will by an appro- 
priate modification of the magnetic field strength. 

Finally, the entanglement properties between the vari- 
ous degrees of freedom of this relativistic system has also 
been accomplished. This hybrid system consists of a cou- 
ple of continuous-variable degrees of freedom associated 
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to the chiral modes, and a discrete-variable associated 
to the spin, therefore constituting a tripartite system. 
We have shown that an entanglement measure, the von 
Neumann entropy of the reduced states, over all possi- 
ble bipartitions can be obtained for the system ground 
state. At the critical point, the quantum correlations of 
the chiral degrees of freedom diverge, which can be also 
considered as a sign of quantum phase transitions. 
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APPENDIX A: SU(1,1) COHERENT STATES 

In this appendix we review the properties of two-mode 
SU(1,1) coherent states. Such states arose from a group 
theoretical approach to the generalization of the usual 
Glaubler coherent states |3ll. which can be constructed 
for an arbitrary Lie group [33 • In this case, we shall 
restrict our attention to the SU(1,1) group, whose gener- 
ators {Kq, K+, have the following Lie Algebra 

[Xo,^±]=±^±, [K_,K+]^2Ko. (Al) 

The Schwinger representation of is algebra (|Aip con- 
sists on a pair of commuting bosonic harmonic oscilla- 
tor operators ar,a; related to the group generators as 

follows [H,!!! 

Ko = \ {a\ar + a\ai + l) , A'+ = a\a\ = (A'_)t. (A2) 

This two-mode representation is very useful, since one 
can use the disentangling theorem [S^] to calculate the 
SU(1,1) coherent states \z,ni) = [/^|ni)|vac)r in Eq. ([20|) 
in the usual Fock basis \n^^ni) 

Jj\ _ g2(-?f+-Ar_) _ gtanhz_R'+glog(cosh"^z)Ko/2g-tanli2i<'_ 

(A3) 

Such factorization of the exponential (|A3[) allows the fol- 
lowing expression for the SU(1,1) coherent states 



\z,ni)=M^, V \ r , ^• (^irtanh'"|z||m + n;,m), 

(A4) 

which has been used in Eq. ([20]) of the left-handed regime 
^ > ^. Analogously, one can obtain the corresponding 



SU(1,1) coherent states in the right-handed regime, using 
a similar disentangling theorem 



\z.n,) =Mn^ Y \ "^• (-l)"tanh"|g||m+n„m), 

(A5) 

which was used in Eq. 

In the Lie algebraic language, the Casimir operator C 
commutes with the generators {Kq,K±} of the algebra. 
In the SU(1,1) case, the Casimir operator turns out to 
be the z-component of the orbital angular momentum 
C — Lz, which is therefore left unchanged under the 
action of the squeezing operator [U^,Lz] = 0. This fact 
has been used in section |Vl where the expected value of 
such Casimir operator is used as an order parameter that 
reveals the existence of a quantum phase transition. 

SU(1,1) coherent states present several remarkable 
non-classical properties, which attain an special interest 
in Quantum Optics since they can be engineered in ex- 
periments by means of a degenerate parametric ampli- 
fier. In the relativistic setting, they arise naturally as 
the eigenstates of a relativistic oscillator subjected to a 
magnetic field along the z-axis. In particular, we shall 
be interested in the following properties of such states, 
which have been extensively studied in [33| . 

{nr)\z,n,) = (ni + l)sh^|z|, 
{ni)\z,ni) = ?^^ch^|^| +sh^|z|, 

{ni)\z,n,) = im + lfsh'^\z\ + (1 + n^ch^lzlsh^lzl, 

{nf)iz,n,) = n\c\iM + sh^lzl + (1 + 3ni)eh2|z|sh2|z|, 

(A6) 

for the SU(1,1) coherent states in the left-handed regime, 
where sh|z| :— sinh|z|,and ch|z| :~ coshjzl, and we use 
the chiral number operators := aJ^Or-, and ni := a|a; 
. Analogously, the corresponding coherent states in the 
right-handed regime fulfill the following 

(^i)|2,n,) = (1 + nr)sh^|z|, 

(^r)|£.n,) rirCh^l^l + sll^|z|, 

(^?>|i,ft.> = (^r + l)^sh'^|i| + (1 + hr)c\?\~z\s\^\~zl 

{nl)\i n ) n^ch^'lzl -I- sh-^l^l + (1 + 3n^)ch2|z|sh^|z|, 

(A7) 

where sh|z| := sinh|z|,and chf := cosh|z|, and we use 
the chiral number operators fij. := dldr, and fii := djdi. 
These expressions (|A6p . (|A7|) are extremely useful for the 
study of the relativistic phonon statistics in Section IVl 



[1] S. Sachdev, Quantum Phase Transitions, (Cambridge [2] N. Lambert, C. Emary, and T. Brandes, Phys. Rev. Lett. 
University Press, Cambridge,England, 2000). 92, 073602 (2004). 



13 



[3] C. Emary, and T. Brandes, Phys. Rev. Lett. 90, 044101 [18 
(2003). 

C. C. Gerry, Phys. Rev. A 37, 3619 (1988). [19 

C. C. Gerry, J. Kiefer, Phys. Rev. A 41, 27 (1990). 
A. P. Mines, C. M. Dawson, R. H. McKenzie, and G. J. [20 
Milburn, Phys. Rev. A 70, 022303 (2004). 

M. Moshinsky and A. Szczepaniak, J. Phys. A 22, L817 [21 
(1989). 

D. Ito, K. Mori, and E. Carrieri, Nuovo Cimento 51 A, 
1119 (1967). [22 
W. Greiner, Relativistic Quantum Mechanics: Wave 
Equations, (Springer, Berlin, 2000). [23 
A. Bermudez, M. A. Martin-Delgado, and E. Solano, 
Phys. Rev. A 76, 041801(R) (2007). [24 
A. Bermudez, M. A. Martin-Delgado, and E. Solano, 
Phys. Rev. Lett. 99, 123602 (2007). [25; 
R. Banerjee, and S. Gosh, J. Phys. A: Math. Gen. 31, [26 
L603 (1998). 

M. M. Horta Barreira, and C. Wotzasek, Phys. Rev. D [27 
45, 1410 (1992). 

G. V. Dunne, R. Jackiw, and C. A. Trugenberger, Phys. [28 
Rev. D 41, 661 (1990). 

E. T. Jaynes and F. W. Cummings, Proc. IEEE 51, 89 [29 
(1963). [30 
D. Leibfried, R. Blatt, C. Monroe and D. Wineland, Rev. 
Mod. Phys. 75, 281 (2003). [31 
Note that these creation-annihilation operators do not 
represent normal modes since they do not commute with [32 
each other, but rather fulfill [ai,aj] = Sijfi, [ai,at] = [33 

<5,,/i, where ^ =: i f ^ - |) , and ^ = i f # + |^ [34 



S. M. Barnett, and P. L. Knight, J. Opt. Soc. Am. B 2, 
467 (1985). 

A. Osterloh, L. Amico, G. Falci, and R. Fazio, Nature 
416, 608 (2002). 

G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev, Phys. 
Rev. Lett. 90, 227902 (2003). 

M. Nielsen, and I. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

A. Galindo, and M. A. Martin-Delgado, Rev.Mod.Phys. 
74, 347 (2002). 

C. H. Bennet, D. P. Divicenzo, J. A. Smolin, and W. K. 
Wooters, Phys. Rev. A 54, 3824 (1996). 

H. Casini, C. D. Fosco, and M. Huerta, J. Stat. Mech. 
0507, P007 (2005). 

P. Calabrese, and J. Cardy, J. Stat. Mech. P06002 (2004). 
L. Amico, R. Fazio, A. Osterloh, and V. Vedral, 
'arXiv:quant-ph/0703044' (2007). 

M. B. Plenio, and S. Virmani, Quant. Inf. Comp. 7, 1 
(2007). 

J. Eisert, and M. B. Plenio, Int. J. Quant. Inf. 1, 479 
(2003). 

G. S. Argawal, Phys. Rev. A 3, 828 (1971). 

A. Serafini, F. lUuminati, and S. De Sienna, J. Phys. B: 

At. Mol. Opt. Phys. 37, L21 (2004). 

R. J. Glauber, Phys. Rev 130, 2529 (1963); 131, 2766 

(1963). 

A. M. Perelomov, Commun. math. Phys. 26, 222 (1972). 

C. C. Gerry, J. Opt. Soc. Am. B 8, 685 (1991). 

L. Gilles, and P. L. Knight, J. Mod. Opt. 39, 1411 (1992). 



